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ABSTRACT 

We study the galaxy environment of active galaxies, radio-loud and radio-quiet quasars 
in the redshift range 0.1 < z < 0.25. We use APM galaxies in order to explore the 
local galaxy overdensity and the bj — R colour distribution of neighbouring galaxies of 
these target samples. For comparison, we perform similar analysis on samples of Abell 
clusters with X-ray emission, and samples of Abell clusters with richness R — 1 and 
R = 0. The projected cross-correlations show that the samples of quasars and active 
galaxies reside in regions of galaxy density enhancements lower than those typical of 
R = clusters. We also find that in the nearby universe the local galaxy overdensity of 
radio-loud and radio-quiet quasars are comparable. The analysis of the distribution of 
bj — R galaxy colour indexes suggests that the environment of quasars is not strongly 
dominated by a population of red galaxies, characteristic of rich Abell cluster , an 
effect that is more clearly appreciated for our sample of radio-loud quasars. 
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1 INTRODUCTION 

Studies of the environment of quasars have been the subject 
of several works in recent years. Some of the analysis in 
the local universe, z < 0.6, suggest that quasars reside in 
groups (Fisher et al 1996), or in clusters (Smith et al 1999, 
Mc. Lure & Dunlop 2000). However, Martinez et al. (1999) 
using a quasar-galaxy cross-correlation analysis suggest that 
high density environments may not be representative of the 
typical quasar neighbourhood at low redshifts. 

At higher redshifts, z > 0.6, some authors have found 
evidence of a different environment for radio-loud and radio- 
quiet quasars (Yee & Green 1987, Yee 1990). Radio- loud 
quasars being associated mainly with rich Abell clusters 
while radio-quiet quasars residing either in groups or in the 
outskirts of clusters. However, more recent studies by Wold 
et al. (2000, 2001) and Best (2000) cast considerable doubts 
on previous results since these new analysis do not show any 
appreciable difference between the typical galaxy overden- 
sity around radio-loud and radio-quiet quasars. 

The discussion about the influence of quasar activity on 
the formation and evolution of galaxies has continued during 
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the last decades. The possibility that the surrounding galax- 
ies could be responsible for triggering and fueling the quasars 
was first suggested by Toomre & Toomre (1972). Interacting 
galaxies and mergers are efficient means of transporting gas 
into the inner regions of the host galaxy. On the other hand, 
the large range of influence of a quasar may have implica- 
tions for structure formation. According to Rees (1988) and 
Babul & White (1991) quasars may ionize the surrounding 
medium up to several megaparsecs away. A different ap- 
proach is adopted by other authors (see for instance Voit 
1996 and Silk & Rees 1998), where a fraction of the energy 
released by quasars is transferred via energetic outflows to 
the gaseous medium. 

In this paper we attempt a characterization of the en- 
vironment of quasars and active galaxies in the nearby uni- 
verse by considering the galaxy density enhancement using 
cross-correlation analysis with APM galaxies. We also use 
the distribution of bj — R colours of neighbouring galax- 
ies to explore their stellar population properties and deepen 
our understanding of nuclear activity and galaxy formation 
and evolution. In section 2 we describe the data, section 3.1 
deals with the cross-correlation analysis, section 3.2 with 
the relative distribution of bj — R colour indexes of APM 
galaxies around the different target samples. Finally, section 
4 provides a brief discussion of the main results. 
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Figure 1. Redshift distributions of the different samples. The 
smooth curve represents the expected redshift distribution for 
APM galaxies with limiting apparent mag nitude b 1 ™ = 20.5 a.s 
derived by Baugh & Efstathiou (1993). 



2 DATA 

The active object target samples used in this work consist 
in radio-loud quasars, hereafter RLQs, radio-quiet quasars, 
hereafter RQQs, and active galaxies, hereafter AGNs, from 
the Quasars and Active Galactic Nuclei 9nth Ed. Catalog 
(Veron-Cetty & Veron 2000). The cluster targets consist 
in clusters of galaxies taken from the X-ray Bright Abell 
Cluster Sample (XBACS) (Ebeling et al. 1996) and Abell 
clusters with richness R = 1 and R — selected from Abell 
et al. (1989). 

All target samples have been restricted to the red- 
shift range 0.1 < z < 0.25 taking into account the 
similarity of the target redshift distributions within this 
range, and the large amplitude of correlation of the tar- 
gets with APM galaxies. Tracer galaxies within 0.5° of 
the targets were taken from the APM South ern Sky 



Catalogues (http://www.ast.cam.ac.uk/~apmcat/) and re- 
stricted to those targets for which data in the two bands bj 
and R are available. The final number of targets in each sam- 
ple are listed in Table 1. The samples of Abell clusters were 
restricted to richness R < 1 to avoid the inclusion of clus- 
ters already considered in the XBACS sample since XBACS 
include mostly R > 2 clusters. 

We have also considered 30 randomly selected fields to 
make an appropriate comparison with the mean background 
properties. Main characteristics of the target samples anal- 
ysed are also listed in Table 1. 

The different target samples have similar redshift dis- 
tributions as it can be seen in Figure 1, so that the results of 
our analysis are not subject to different 7f-correction effects 
as well as other possible systematic that depend on sample 
depth. For comparison, in each panel of Figure 1 is shown 



Figure 2. Projected Cross-Correlation Functions. Filled circles 
correspond to XBACS cluster sample, filled squares correspond 
to Abell R = 1 and filled triangles to Abell R = 0. The best fit 
power-law to the projected cross-correlation function for the clus- 
ter samples are shown (see parameters in text). Empty squares 
correspond to RLQs, empty triangles to RQQs and empty circles 
to AGNs. 



the expected redshift distribution of APM tracer galaxies 
for a limiting apparent magnitude 6j m = 20.5 (see Baugh & 
Efstathiou 1993) in arbitrary units. 



3 ANALYSIS 

3.1 Galaxy overdensities around the targets 

In order to make a suitable characterisation of the typical 
galaxy density enhancement around the objects in the se- 
lected samples, we computed the cross-correlation functions 
for centers of the different target samples and APM trac- 
ers with a limiting magnitude 6j m = 20.5. Since all center 
targets in our samples have measured redshifts, we are able 
to compute cross-correlations as a function of projected dis- 
tance which improves the information obtained from angular 
correlations. 

We have used the following estimator of the two-point 
correlation function: 



W(r P ) = 



DD(r p )N R 
DR(r v )N D 



(1) 



where DD(r p ) and DR(r p ) are the number of target-galaxy 
and target-random point pairs respectively, separated in 
each case by a projected distance r p in the interval [r p — 
Ar p ,r p + Ar p ], Nn and Nr are the number of objects in 
the data and random catalogs respectively. When comput- 
ing projected distances we have assumed a flat cosmology 
(go = 0.5) and a Hubble constant Hq = 100/ikms _1 Mpc _1 . 
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Table 1. Target sample characteristics. 



Sample 


Number of objects 


Source 


Main characteristics 


RQQs 


41 


Veron-Cetty & Veron (2000) 


Mb < —23, radio-quiet 


RLQs 


19 


Veron-Cetty & Veron (2000) 


Mg < — 23, radio-loud 


AGNs 


35 


Veron-Cetty & Veron (2000) 


M B > -23 


XBACS 


34 


Ebeling et al. (1996) 


Fiim = 5.0 X 10~ 12 erg cm -2 s _1 


Abell R = 1 


35 


Abell et al. (1989) 


R= 1 


Abell R = 


45 


Abell et al. (1989) 


R = 



The resulting cross-correlation functions are shown in Fig- 
ure 2. Error bars in these figures were estimated with the 
bootstrap resampling technique (Barrow et al. 1984) with 
30 bootstrap samples. 

In Figure 2 we show the projected cross-correlation 
functions corresponding to the three cluster samples and 
the samples of active objects, where it can be seen the de- 
crease of the cross-correlation amplitude moving from richer 
to poorer clusters as well as the lack of strong correlation for 
the samples of active objects. It is clear that the projected 
cross-correlation functions of the samples of active objects 
have systematically lower amplitude than those of cluster 
samples (even R = 0) and they are indistinguishable taking 
into account uncertainties. 

Assuming a power-law model for the spatial correlation 
function, £(r) = (r/ro) the projected correlation function 



Ntarget = (N) — Nran 



(•») 



W(r p ) = Ctt 



1/2 r(( 7 -l)/2) 



r(7/2) 



(7- 



where the constant C, 



C 



(2) 



(3) 



takes into account the probability p(x) for a tracer galaxy to 
be found at a distance x from the observer and is derivable 
from the galaxy luminosity function 4>(L): 



p(x) oc 



<j>{L)dL, 



(4) 



where a galaxy at a distance x must have a luminosity 
greater than L m i n (x) to be included in the tracer sample 
(see Lilje & Efstathiou 1988). In equation (^) yi is the ra- 
dial distance to the i—th center target. We have used a 
Schechter fit to the luminosity function (Schechter 1976) for 
the tracer galaxies with parameters M£ = —19.7 + 5 log ft, 
q = —1.0 and we have included a linear A"-correction, 
K = 3z, in computing L m i n (a;). By fitting power-laws to 
the cross-correlation functions we have inferred slopes and 
correlation lengths for the cluster samples. We obtain 7 = 
2.09±0.09, r = (ll±2)/T 1 Mpc for XBACS; 7 = 1.9±0.2, 
r = (11 ± 3)/i _1 Mpc for Abell R = 1; 7 = 1.9 ± 0.2, 
r = (9 ± 2)/i" 1 Mpc for Abell R = 0. 

The low amplitude and the large error bars of the cross- 
correlations of active objects do not allow adequate power- 
law fits. A useful quantity to characterise the typical galaxy 
overdensity around the center targets considers the number 
of tracer galaxies at projected distances r p < l/i _1 Mpc with 
respect to the random fields, 



Our sample of tracer galaxies are limited to b.j < 20.5 and 
corresponds to objects brighter than Mbj = — 17.7 + 5 log h 
at z ~ 0.15. We find for active objects NXgn 3 ~ 15; 
NrqqI ~ 13; N c r Tq! ~ 10. For clusters we find N^bTcs ~ 
100; ATn=r ~ 50 and A^g ss ~ 30. From these results it is 
clear that active objects are not typically located in moder- 
ately rich galaxy enhancements such as Abell R = clusters, 
as suggested in previous works (see for instance Smith et al. 
1999, Mc. Lure & Dunlop 2000). 



3.2 Distribution of b.j — R colour indexes of 
neighbouring galaxies 

The relative fraction of different stellar populations in galax- 
ies are reflected in the bj — R colour index. Elliptical galaxies, 
dominated by an old population of stars, have typical values 
bj — R ~ 1.5 — 2, while in late types with recent episodes 
of star formation, b.j — R < 1. In this subsection we explore 
the properties of the stellar populations of galaxies in the 
vicinity of the different target samples using b.j — R colour 
indexes. We study the normalised distribution of bj — R of 
APM galaxies around the targets within a projected distance 
r p < 0.5/i _1 Mpc. We adopt the same limiting apparent mag- 
nitude bj m = 20.5 for all tracer samples and we compute: 



f(h ~R) = 



N(bj - R) 
N TOT A{bj - R)' 



(6) 



where N(bj — R) is the number of APM galaxies with 
colour indexes in the interval defined by [(bj — R) — A(bj — 
R)/2,(bj-R)+A(bj-R)/2], and N TOT is the total num- 
ber of APM galaxies within rp < 0.5h~ Mpc. Error bars 
in these distributions were computed using the bootstrap 
resampling technique. We subtract from these distributions 
the colour distribution of tracer galaxies in the randomly 
selected fields. The results are displayed in Figure 3. As it is 
expected from the well known morphology-density relation 
(Dressier 1980), there is a marked relative excess of early 
(1 < bj - R < 2) to late type galaxies (0 < bj - R < 1) 
within r p < 0.5/i _1 Mpc for the three samples of clusters, an 
effect that is stronger for the richer cluster samples. Active 
objects show a weaker tendency, although, the distribution 
of colours of RQQs shows a similar excess of red galaxies 
than in R = clusters. For RLQs, the signal for the pres- 
ence of a red galaxy population is almost negligible, while 
AGNs show an intermediate behavior between RQQs and 
RLQs. 

To quantify the previous results, we have computed a 
parameter V that characterizes the relative fraction of red 
to blue galaxies in a given range of projected separations 
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Table 2. V parameter for the different target samples within three projected distances. 



Sample 


V 

r p < 0.3ft -1 Mpc 


V 

r p < O^h^Mpc 


V 

r p < l/i _1 Mpc 


RQQs 


0.10 ±0.06 


0.09 ± 0.04 


0.07 ±0.02 


RLQs 


-0.22 ±0.08 


-0.05 ±0.07 


-0.05 ±0.04 


AGNs 


0.08 ± 0.05 


0.03 ± 0.04 


-0.01 ±0.02 


XBACS 


0.49 ± 0.03 


0.38 ± 0.02 


0.30 ±0.02 


Abcll R = 1 


0.17 ±0.05 


0.20 ±0.04 


0.16 ±0.02 


Abell R = 


0.18 ± 0.07 


0.09 ± 0.04 


0.06 ±0.02 


Field 


-0.11 ±0.01 
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Figure 3. Colour index bj — R distributions normalised to 
randomly selected fields corresponding to galaxies within r p < 
0.5h~ 1 Mpc from the targets. 



N R + N b , 



(7) 



where Nr is the number of neighbour tracer galaxies in 
the range 1 < bj — R < 2 and Nb, is the number in the 
range < bj — R < 1. In Table 2 we list this parame- 
ter for three projected distance ranges (r p < 0.3/i _1 Mpc, 
r p < 0.5/i _1 Mpc and r p < l/t -1 Mpc) and their correspond- 
ing bootstrap resampling uncertainties estimates. 

The values of V in Table 2 allow us to quantify and 
extend the previous results of this subsection. By inspection 
to this table it can be seen that galaxies in clusters show 
the tendency of systematically larger colour indexes near 
the center. The population of galaxies in the neighbourhood 
of RQQs have V comparable to that of galaxies in R — 1 and 
R — Abell clusters. It is worth to notice the negative values 
of V for galaxies in the vicinity of RLQs which indicate the 
lack of a red galaxy population, characteristic of high density 
regions, around these objects. In order test the stability of 
these results we have changed the boundaries adopted to 



calculate the values of V ■ The analysis shows that the above 
conclusions are robust and do not depend on the particular 
choice of colour boudaries in the calculations, provided the 
parameter V reflects the relative fraction of early-type to 
late-type galaxies. 



4 DISCUSSION 

The results discussed previously lead us to conclude that 
quasars and AGNs are not placed in galaxy density enhance- 
ments corresponding to rich and moderately rich clusters of 
galaxies. As it can be clearly appreciated by inspection to 
Figure 2, the galaxy overdensity around these center tar- 
gets is significantly lower than that corresponding to Abell 
R — 1 or R — clusters of galaxies. Thus, active objects 
are more likely to reside in groups (typically 10 neighbours 
with M bj < Af 6 * 7 + 2 within r p < l/i _1 Mpc). It has been 
often argued that quasar environment correspond to that 
of moderately rich clusters, our results indicate that active 
nuclei reside in regions of at most half the projected galaxy 
density of R = clusters. 

The distribution of bj — R colour indexes of galaxies 
within projected distances r p < l/i _1 Mpc centered on RLQs 
is significantly different than that corresponding to clusters. 
This fact suggests that the stellar population of the neigh- 
bouring galaxies of RLQs differs significantly from that in 
rich and moderate galaxy density enhancements. The envi- 
ronments of RQQs show bj — R colour indexes similar to 
galaxies in clusters while AGNs show an intermediate be- 
havior. 

We find statistical evidence that the distribution of b j — 
R colours of galaxies in the vicinity of radio-quiet and radio- 
loud quasars are different. The latter have a population of 
neighbouring galaxies with lower values of colour indexes. 
This fact should be considered in models of the effects of 
nuclear activity on galaxy formation and evolution. 
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